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The absence of fermionic, asymptotical one-particle states in the Luttinger model raises the sus- 
' picion that the interactions are actually strong at the vicinity of the Fermi points. The functional 

I internal space renormalization group method, a systematical scheme for the computation of effective 

, coupling strengths off the Fermi points, is applied to shed some light on this issue. A simple trunca- 

tion of the evolution equation shows that the theory is indeed strongly coupled at the Fermi points 
for arbitrarily small value of the bare coupling constant and the removal of the cutoff is blocked 
, , by Landau poles. A peculiar feature of the normal ordering scheme is pointed out to explain the 

■ absence of these effects in the bosonized solution. 
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I. INTRODUCTION 



One dimensional massless fermions display rather peculiar features. Owing to another specialty of one dimension, to 
^ the equivalence of the helicity and the chirality, the left and right moving fermion numbers are individually conserved 
in chiral invariant theories, like in the massless Thirring [1] or Gross-Neveu [2] models. The physics of electrons 
moving in a one-dimensional wire can effectively be described by linearizing the dispersion relation around the Fermi 
points. The resulting Luttinger model [3], a theory of relativistic, massless fermions, is the massless Thirring model 
] except that the dispersion relations of the left and right moving fermions are displaced horizontally with respect to 
(~| . each other in order to accommodate the chemical potential. This model provides a rich structure to compare with 
O ' experimentally well controlled situations [4]. 

^ O^ , One dimensional models can elegantly be renormalized by means of normal ordering and the works aiming at 
the Luttinger model are almost exclusively based on this method. Most widespread use of normal ordering is the 
T-H \ bosonization. The bosonization of the excitations within a given fermion number sector of the Fock space is in 
J> ■ principle possible in any dimensions. But, another interesting low dimensional "accident", the bosonized effective 
' theory [5-7] is local in dimension. The bosonized Luttinger model consists of free bosonic excitations which stand for 
2^ ' the asymptotical, multi particle-hole states of the fermionic representation [5, 6, 8]. The most striking phenomenon of 
' this chiral model is the loss of Fermi-liquid properties, the continuity of momentum distributions at the Fermi points 
. at vanishing temperature. In the absence of a transverse direction to the motion the slightest interactions destabilize 
' the particle-like asymptotical states which always propagate with the same speed. The asymptotical states of the 

theory, if they exist, must correspond either to massive particles or to non-interacting quasi-particles. 
' The non- interacting nature of the asymptotical sector seems to be consistent with the vanishing of the beta functions, 
' confirmed in a number ways, independently of the bosonization. One chain of arguments is based on the application 
^ ■ of Ward identities which are used to show a perfect cancellation between the self energy and the vertex corrections 
and to render the simple RPA exact [9-11]. The Wilsonian blocking employed in momentum space for the coupling 
constant of the bare action produces vanishing beta function [12, 13], too. 

It is worthwhile mentioning a superficial analogy with the confinement of color charges in QCD. In fact, the non- 
^ , propagating nature of the elementary fermionic excitations in one dimension is equivalent to the cancellation of the 
i;^ ■ residue of the particle pole in the Green functions and the scattering amplitude. These are in turn the hallmark of 
confinement in hadron physics, the absence of colored asymptotical states. But the closer look betrays that the two 
phenomena are actually different. In fact, the bare particles of the one-dimensional world remain non- propagating for 
^ ] either sign of the coupling constant, i.e. both for attractive and for repulsive interactions, as a result of the kinematical 
origin of the effect. On the contrary, the attractive color forces are crucial in removing the colored asymptotical states 
in QCD. 

Nevertheless, the superficial analogy between these cases raises other puzzles. The quark confinement is already 
a rather spectacular phenomenon because it changes the degrees of freedom with the scale of observation. But the 
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strong correlations between far separated quarks is reflected in the strong interactions among the quasiparticles, the 
hadrons. What happens in the bosonized model is even more surprising: How could the complete screening between 
"hadrons" take place in contrast to the strong, long-range correlations between the elementary constituents? Another 
question comes from the fact that classically scale invariant, renormalizable models display dimensional transmutation 
[14] . It arises when a (partial)resummation of the perturbation expansion produces a combination of the cutoff A and 
the dimensionless bare couplings in observables which remains finite in the limit A ^ oo. This must happen unless the 
IR scaling laws remain the same as in the UV regime. The non-interacting, massless particle-hole quasiparticles posses 
no scale, any scale generated would indicate finite interaction strength. Why there is no dimensional transmutation, 
or different scaling regimes in the Luttinger model? No other classically scale invariant, renormalizable model shows 
such a simple structure. 

These questions lead to a one-loop renormalization group study of the Luttinger model [15] which gave non-trivial 

beta functions, strongly coupled dynamics and singular momentum dependence in the effective interaction strength 
around the Fermi point. Furthermore, there were UV Landau poles in certain kinematical regimes of the effective 
coupling strengths, suggesting problems with renormalizability. This introduces the question whether the Luttinger 
model is indeed renormalizable. 

We argue below that renormalization schemes based on normal ordering are rather fragile in the sense that cutoff 
independent, finite physical effects, suc;li as dimensional transmutation and Landau polos, may arise from parts of the 
renormalized action which are vanishing when the cutoff is removed. In order to see if such phenomena are indeed 
present in the Luttinger model we turn to another powerful method, the functional renormalization group in the 
internal space [16]. 

In the usual renormalization group strategy the UV cutoff is lowered and the dynamics of the degrees of freedom, 
eliminated in this manner, is taken into account by adjusting the parameters of an effective bare action. On contrary, 
the fluctuations are turned on gradually in the internal space scheme and their effects arc accumulated by the 
evolution equation for the effective action [17]. The advantage of this method is that it induces weaker non-local 
artifacts as blocking in momentum space and it gives access to effective interaction strengths off the Fermi points. 
These features are essential for the Luttinger model where eventual strong interactions appear as a singularity of the 
effective interaction strengths in momentum space in the vicinity of the Fermi points. 

The only approximation needed in the functional renormalization group method is the restriction of the evolution 
equation onto a functional space with manageable size. The perturbation expansion and the functional renormalization 
group based solutions of a model can be contrasted by stating that one resums infinitely many orders generated by 
few relevant vertices in the former and uses one-loop expressions but infinitely many effective vertices in the latter 
procedure. We made a modest step in this work in extending the computation towards more coupling constants and 
considered the Luttinger model with two running coupling constants whose flow is obtained in the one-loop level. 

Our results conflrm the absence of evolution of the coupling constant when it is considered just at the Fermi points. 
But as soon as one comes off the Fermi points even by an infinitesimal amount our renormalization group scheme 
produces a non-trivial evolution which renders the theory strongly coupled in the vicinity of the Fermi point for 
arbitrary bare parameters. The cutoff generates Landau poles which seems to block the way of the removal of the 
cutoff. These results contradict the picture emerging from the normal ordering solutions and suggest the survival of 
cutoff effects. 

The plan of the paper is the following. The model and its generating functional are introduced in Section II. A 
peculiar feature of the normal ordering schemes is discussed in Section III. The internal space renormalization group 
is briefiy described and applied to the Luttinger model in Section IV. The projection of the evolution of the equation 
into a suitable chosen functional subspace is discussed in Section V. The projected evolution equation for the effective 
action is given in Section VI. The simplest application of the evolution equation is presented in Section VII where the 
effective coupling strength is chosen at the Fermi points and the vanishing beta function is recovered. Section VIII 
contains the main results of this work, the integration of the beta functions for coupling strength off the Fermi points. 
Our results are summarized in Section IX. Two appendices facilitate the reading of the paper. The first sums up the 
notation and conventions used in the paper and the second contains few details of the computation of the one-loop 
integrals of the evolution equation. 



3 



II. LUTTINGER MODEL 



The Luttinger model and its generating functional for the connected and the one-particle irreducible (IPI) Green 
functions are introduced in this section. The model is defined by the action Sltjj^ , ip] = S'o[V'^ V'] + Si[ip^, tp], with 

(Ti,<T2,0-3,0-4 "'2:i,1/l,X2,y2 

where the field operator ip^-^ stands for left or right moving fermions with cr = — 1 or +1, respectively and 

e<7(p) =<jp-kF. (2) 

Sec Appendix A for notations and conventions. The separate conservation of number of left and right movers is 
achieved by allowing the interactions gtl%'^xt,v2^ 9x;^'~x2,y2^ 9t;tiuX2.y2^ ^"^^ 9x;^[txt,y2 °^^y ^° non- vanishing. 
Symmetry under space inversion requires g+;+[+^+y^ = 9~lyQ^^y^ and 9tityui2,v2 = 9xl~yuxt,y2- 
In order to determine the space-dependence of the coupling constants one introduces the Fourier transforms 



and 

Cri,Cr2,<T3,0-4 _ / „ipiXl-ip2yi+ip3X2-iy2Pi X ;^CTl,CT2,0-3,Cr4 (A\ 

yXi,Vx,X2,y2 ~ ^ "Pl-P2+P3-P4,;0ypi,p2,p3,p4 ' 

•'Pl,P2 ,P3 ,P4 

The coupling constants which order the Grassmanian fields in the path integral according to Eq. (1) are given by 

= iri(2p°-pl-pl+2p°) X 

^Pl,P2,P3,P4, y 

= ir,{2p°-p°-p°+2p'i) n /gx 



where and g" are constants and 77 = 0+. 

The generating functional for the connected Green functions is given by 

The effective action, defined by the Legendre transformation, 

r[V;+, V] = -W\j,3^] + E • + jl ■ (7) 

(T 

with 

(t;x 



5W 



= t:^ (8) 



generates the IPI vertex functions by its functional Taylor expansion. 

The generating functional (6) needs regulator. For the sake of a better insight into the cutoff effects let us mention 
four easily implementable regulators, (i) The simplest regulator is a sharp cutoff in momentum space, the restriction 
of the Fourier modes of the field configuration to |p — crfci?| < A. (ii) We can construct a slightly different sharp cutoff 
regulators by placing the system onto space-time lattice. The most important change is the bending of the dispersion 
relations into smooth, periodic functions, (iii) We shall need smooth cutoff which is easiest to reach by smearing the 
local terms in the interaction Lagrangian. One possibility is to introduce the smeared field 



V'B a,x = / Px-yi^a,v (9) 

Jy 



4 



and use the action Sltp^, tp] = 5o[V'^, '^i'] + 5'i[V'B) V's] i^)- ^ natural choice for smearing is 

P. = -e-^'^' (10) 

TT 

which amounts to the use of 

V'B <r,p = Pp'4}a,p (11) 

with 

Pp = e-^'/'^\ (12) 
We note that the smearing function pp — pa{p)j with 

brings us back to sharp momentum cutoff, (iv) Another smooth cutoff which preserves the local form of the current 
in the interactions and sphts the two currents only is obtained by using the interaction 

in Eq. (6). Each regularization admits bare Ward-identities. 



III. BPHZ SCHEME AND BOSONIZATION 



Most of the studies of the Luttingcr model was carried out in an elegant subtraction scheme in the Hamiltonian 
formalism which is based on normal ordering. We highlight in this section the peculiarity of the results by showing 
that the normal ordering is actually very similar to the BPHZ subtraction scheme. 

The BPHZ subtraction scheme is based on the observation [18, 19] that the primitive divergences of an irreducible 
Green function appear in the first few order of the Taylor expansion in the external momenta around zero. These 
divergences are subtracted by the introduction of local counterterms, given just by these few divergent terms. This 
procedure can be organized in such a manner that the subtraction is carried out on the level of the integrand of the 
loop integrals [20]. The resulting scheme is remarkable because it has no divergences whatsoever, the loop integrals 
are manifestly finite [21]. 

The bosonization in 1+1 dimensions has been thoroughly discussed in a number of papers [3, 5-7] and we are content 
to make few remarks about the cutoff effects by making the regularization explicit. The bosonization starts with the 
explicit construction of the correspondence between the fcrmionic and the bosonic field operators and continues with 
the demonstartion that the Hamiltonian of the Luttinger model is of massless, noninteracting boson fields [3, 5]. We 
shall consider the the verification of the bosonic commutators in details. 

The field operator corresponding to a chiral fermion in a box of size L and periodic boundary conditions is given 

by 

k 

k = 2imk/L, where Uk is integer. These operators obey the set of canonical anticommutation relation whose only 
non- vanishing member in momentum space is {ck, c^./} = 5k,k'- The boson operators are introduced by the equation 

\ = — 1=Y. 4-9Cfc (16) 

for g > and 

9>0 
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The canonical bosonic commutators for bg and 6j follow in a trivial manner except the relation for 

l^^' = ;r E ci_^ck-g. (18) 



The diagonal matrix elements of this equation contains the difference of two divergent quantities and needs a cutoff. 
The regulator plays two roles simultaneously: It (i) makes the expressions finite by suppressing the short distance 
modes and, and (ii) defines an order of the summation/integration. The latter is rarely needed because the coun- 
terterms, pulled into the loop summation/integration render the convergence absolute and hence independent of the 
order. But there might observables with sums/integrals which are finite but non-absolutely convergent. These have 
no counterterms but still need a definite and consistent prescription in order to render the theory well defined. These 
observables arc called anomalous. 

In lattice regularization the Brillouin zone is periodic and the shift k ^ k + q in the first term is allowed and renders 
the commutator vanishing. The same argument applies for sharp momentum cutoff, as well [3]. 

But the regularization scheme where smooth cutoff is combined by normal ordering gives different result [5]. We 
shall use the regulator (iii) introduced in Section II which yields 

bsq = ^ E P*k-qPkci_gCk (19) 

k 

as bare operators for the interaction Lagrangian. We define the fermionic vacuum by cfe|0) = c|.,|0) = for k, -k' > 
and the corresponding normal ordering of an operator O is 

:O:=O-(0|O|0). (20) 

The regulated Eq. (18) now reads as 

[bBq,big] = ;;J-El^fc-9l^4cfc- :;^El^'=l^4-gCfc-9 (21) 

= ;r E \pk-i\^ 4cfe ---^Yl i/'fci' 4-,cfe-5 : +^ E ipfci'(oi4-,cfe-«io) - E ipk-qfrncWm. 



One can show that any matrix element of the sums is absolute convergent and the reordering k ^ k + q is allowed, 
giving 

[bBq,bl^] = t{l + Cg) + Dg (22) 

where 

^1 k<0 ^1 fe=0 

^Y.^\pk-q?-\pk+q?):clck: (23) 



Dq 



In the works using bosonization the usual regulator is (iv) of Section II which leads to Eq. (21) with = 1. The 
sum with the normal ordered product gives absolute convergent matrix elements but the last two sums diverge and 
renders this cutoff unacceptable. 

In order to understand better the problem of the regularization (iv) let us point out the similarity of the normal 
ordered scheme and the BPHZ method. One subtraction is sufficient in 1+1 dimensions and only its value differs 
in the normal ordered and the BPHZ schemes. The normal ordering shares the advantages of the BPHZ scheme, 
namely simplicity, the disappearance of divergences. But this does not mean that these schemes need no regulator 
whatsoever. The following two circumstances make the regulator unavoidable. 

If there arc non-absolutely convergent but finite loop summations/integration, anomalous observables, in the theory 
then regulator is needed in order to define their value in a consistent manner. The consistency is important, eg. one 
could have chosen any order in the second summation of the last line in Eq. (21) with = 1 which yields finite 
result. But it is not obvious how to define sums which appear in other loops. One needs a consistent regulator which 
settle this issue on the level of the bare Hamiltonian and not of particular expressions. 
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The other instance which requires regulatrization is when the cutoff manages to generate a finite scale. An example 
is the dimensional transmutation. The only way the dimensional transmutation can appear in the BPHZ scheme is 
through the finite cutoff-dependence which is vanishing in the limit A oc because all divergent cutoff dependence 
which generates the new scale in other schemes are now subtracted. Therefore the finite expressions of the cutoff 
which tend to zero in the limit A ^ oo must be retained before resumming the perturbation expansion in order no tot 
miss dimensional transmutation. There are simpler or more complicate schemes, where the same renormalized physics 
is recovered in a straightforward or more involved manner. The BPHZ scheme which is very powerful for massive 
models without dimensional transmutation but becomes highly nontrivial for massless models. This happens because 
its renormalized parameters are not closely related to the true effective coupling constants and finite, physical effects 
(eg. dimensional transmutation) arise from infinitesimal small modifications of the BPHZ renormalized dynamics. As 
far as the BPHZ is concerned, one can save it with the expense of some complications, by expanding the irreducible 
Green functions around non-vanishing energy-momenta, but the normal ordering can not be rescued in the absence 
of any scale to rely upon. 

Another case when the regulator is needed is when the cutoff can not be moved beyond a certain threshold. This 
can easily happen in a perturbatively renormalizable theory, where the divergences are eliminated order by order in 
the perturbation expansion. In fact, the renormalization condition which fixes the theory is given as a set of equations, 
containing certain experimentally observed values on one side and the corresponding expressions on the other side. 
These expressions contain the cutoff and the bare parameters in a highly non-linear manner and it may happen that 
they have no solution with finite values of the bare parameters beyond certain scale, usually called Landau pole scale, 
for the cutoff. 

After these general remark we now return to the discussion of Eq. (22). The commutation relation of Ref. [5] is 
obtained from Eq. (22) by setting Cq = Dq = 0, quantities of 0{q^/A^). After this simplification the bosonization 

is quickly completed by expressing the normal ordered fermionic Hamiltonian by means of the boson field and by 
performing a Bogoliubov transformation. The result is a normal ordered Hamiltonian for a massless, non-interactive 
boson field. 

Our point is that the regulator (iv) regulates all Feynman graphs contributing to Green functions because it 
suppresses interactions beyond the cutoff scale. But it does not remove divergences, such as in Eq. (22), from the 
unperturbed dynamics even if the regulator is amended by normal ordering. Normal ordering is a reordering of the 
loop contributions only and not regulator. The control of all divergences actually leads to the Oicf' /A?) terms in Eq. 
(22). They are certainly negligible as far as the canonical commutation relations are concerned at face value in the 
renormalized theory but they modify the bosonic Hamiltonian which is based on these relations and generate some 
0{h.~^) interactions. These interactions which are dropped by the formal limit A ^ cx) in the normal ordered scheme 
may generate finite effects via dimensional transmutations or Landau poles. 

We shall, in this work, look for cutoff effects in the theory defined by sharp momentum space cutoff in order to 
see if the omission of the (^(g^/A^) terms in Eq. (22) is indeed justified. This choice of regulator brings us back 
to the discussion of different regularizations after Eq. (18). Have we not already lost bosonization with all of its 
consequences when sharp momentum cutoff is used? 

The fact that different, but acceptable regularizations produce different result for a finite value of the cutoff is 
not surprising. All what one expects from renormalizable theories is that different cutoff theories converge to the 
same renormalized one when the cutoffs are removed. The usual argument based on universality assures that there 
is a unique Luttingcr model, as long as its theory is renormalizable, whatever acceptable regulator is used. It might 
naturally happen that certain effects are more or less easily identifiable when a given regulator is used but the 
numerical value of the observables must agree. For example, there is a well known no-go theorem [22] that chiral 
fermions, as in Eq. (15), can not be accommodated within the framework of the usual quantum field theories when 
lattice regulator is used. But rather involved constructions which can be viewed as particular regularization schemes 
and violate the conditions of the no-go theorem can realize chiral lattice fermions in particular limits. In a similar 
manner, we may regard the fragility of bosonization as an indication that this scheme concentrates a large amount 
of informations in the ^^(A"^) part of the renormalized action. In other schemes these informations spread over the 
finite part of the action. If they play appreciable role in the dynamics then bosonized models require an unusually 
careful treatment. 



IV. RENORMALIZATION GROUP IN THE INTERNAL SPACE 

The traditional renormalization group method [23] consists of the successive elimination of the modes. There are 
two possible strategies in determining the order of the elimination. The conventional way is to eliminate modes at the 
running UV cutoff A which moves in the IR direction. One can imagine two different ways to justify this procedure, 
(i) One starts with the unavoidable non-local effects of quantum field theories. It is well known that a theory is 
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always non-local at distances which are shorter than the UV cutoff scale. These involved and dangerous correlations 
are turned on gradually during the blocking if the UV cutoff moves from the UV to the IR direction, (ii) Another 
strategy to fix the order of the modes to be eliminated is that the more perturbative modes should be eliminated by 
means of a simpler effective theory. Hence one should tackle with the more nonperturbative modes by means of a 
more sophisticated, "dressed" effective theory which is obtained after having already eliminated a large number of 
modes. The kinetic energy usually renders the short distance modes perturbative therefore the modes are eliminated 
in decreasing order of their momentum, their scale parameter in the space-time, called external space. This latter 
corresponds to the convention of considering the field variables (/)(a;) as a mapping from the external into the internal 
space. 

The caveat in the reasoning (i) is that the non-local effects usually reach well beyond the naive length scales 1/A. In 
fact, let us consider a smooth momentum space cutoff which leaves the modes with momentum p < A unchanged and 
suppresses modes withp > cA (c > 1). The details of the suppression mechanism appears in effective vertices acting on 
distance scales [1 / cA, 1/ A] . These details of the cutoff arc arbitrary and non-physical therefor these non-local effects 
are non-physical, too. In particular, sharp momentum cutoff always generates infinitely long range correlations. Their 
strength is vanishing as A ^ oo in a renormalizable theory but they are present in effective theories where A is finite. 

One might argue that the subsequent steps of the lowering of the cutoff always cancel the non-local contributions 
of the previous step and replaces them by the non-local contribution corresponding to the lowered value of the cutoff. 
But any truncation of the evolution equation raises the issue of approximate cancellation and forces us to search of 
alternative blocking schemes which generate as weak non-local effects as possible. 

Such an alternative scheme can be found by generalizing the reasoning (ii). The point is that there are two ways 
a mode may become perturbative. One, covered by reasoning (ii), is that the corresponding action is dominated by 
the kinetic term. But another possibility is that its amplitude is suppressed by the potential energy. This second 
way a mode may turn out to be perturbative is exploited when the order of the elimination proceeds according to the 
amplitude of the modes, the scale in the internal space [16]. One starts with modes of small amplitudes leaving the 
large amplitude modes at the end. This procedure amounts to the realization of the Callan-Symanzik proposal [24] 
when the amplitude of the modes are controlled by the momentum independent mass term in a relativistic theory. 

There are cases when the momentum independent suppression is not available, when symmetry excludes mass like 
in gauge theories. A possible solution to this problem can be found by generalizing this scheme and taking any 
continuous parameter in the dynamics as a control parameter which (a) corresponds to a quadratic expression of the 
fields and (b) the theory is soluble for certain limit of this parameter. The condition (a) may require the realization of 
this program within the framework of the effective theory for composite operators, such as (jP in the (f)'^ scalar model 
[16] or fermion bilinears [25-28]. 

What kind of non-local effects are generated by renormalizing in the internal space? This is an essential question 
for the Luttinger model where the collinear IR divergences are expected to plague the perturbation expansion and 
to generate long range effective interactions in terms of the elementary fermionic excitations, as mentioned in the 
Introduction. 

The Callan-Symanzik scheme in which the momentum independent, relativistic mass term is used to control the 
elimination of the modes [16] appears to be a scheme with a minimal amount of non-locality. When the mass term is 
forbidden, as in the Luttinger model, then wc may use the coupling constant as a control parameter. More precisely, we 
rescale the coupling constant by a factor which will be evolved from to 1. This scheme will yield a gradual turning 
on the "dangerous" part of the quantum fluctuations which are not factorizable according to the Wick- theorem. The 
only non-local feature of this suppression scheme comes from the energy- momentum dependence of the interactions, 
needed for the separation of the left and right modes. But this is a A- independent real physical effect and is a necessary 
ingredient. 

Another advantageous feature of this scheme when applied for fermionic systems can be seen in the following 
manner. In fermionic models at finite density the dangerous, soft modes arc around the Fermi surface and the 
blocking should zoom in this area at the final stage of the evolution. This usually requires a rather involved geometric 
construction which can easily be avoided by controlling the interaction strength. A rescaling of the fermionic fields by 
removes the factor A from the usual quartic interaction term and generates the coefficient 1/A in front of the 
quadratic part of the action, i.e. multiplies the propagators by A. One can see again that the momentum dependence 
brought by this control parameter is dictated by the the true, A-indepcndent dynamics, the kinetic energy, instead of 
an artificial device and it may build up the eventually strongly coupled dynamics around the Fermi surface. 

Beyond the issue of preserving locality better, the most important difference between the external and internal space 
schemes from the point of view of fermionic system is the possibility in the latter case to obtain effective interaction 
strengths off the Fermi surface. In fact, the external space blocking methods zoom into the Fermi surface at the end 
point of the evolution only. As a result, any off- Fermi surface region is handled without the dynamical contributions 
of modes between the region in question and the Fermi surface. This happens because these modes appear at a later 
stage of the evolution only when the effective coupling strengths off the Fermi surface have already been fixed. On 
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the contrary, the internal space scheme blocking has no impact on the external space and leaves the possibility of 
tracing the evolution of the effective coupling strengths at in an arbitrary kinematical region. 
The realization of the internal space blocking scheme starts with the extension 



„cri,o-2,o-3,CT4,/,t ,/, J,\ J, 

yXi,Vi,X2,y2 ^'ri,Xi Vcr2,yi r(T3,X2 



(24) 



<7i,<72,<^3,<^4 ■' Xl,yi,X2,V:i 

of the action (1) which can be written after the rescaling ip X^/'^ip, -0^ A^/^-f/'t as 



^ „ Jx „, „. „. J xi,yi,X2,y2 



a\ .(72 -(T3 ,(74 



The evolution equation for the generating functional is obtained by simply applying the derivative with respect to 
the control parameter A on Eq. (6) where the action is given by Eq. (25), 

X l^^l;x[-dt + ea{-id^Ma;x 

= ^ E [(^t + Wjt^i^) ■ [-^t + ^<^i-i9x)]] . (26) 

(7 

The truncation of this functional differential equation is rather difficult because the functional W[ja,jl] is too com- 
plicated due to strong non-local features. The IPI vertex functions are supposed to be more local than the connected 
Green functions therefore we express this equation in terms of the effective action. We introduce the notation 



^(2) ^ 6^F[(^,x] 



(27) 



use the relation 



.(2) 



.(2) 



p(2) 

r(2)' 



.(2) 
-(2)' 



p(2) \ 

+ ,2' +;y ' 



.(2) 
.(2)' 



,-(2) 



-5z,x 



r(2) 

/l 0\ 

10 
10 

Vo 1/ 



.(2) 



p(2) 
p(2)' 
T-(2) 



-wf . 

\ J-,y,J+;x 



3+ y^3- 



w. 



(2) 



3+,y,3- 
(2) 



(21 



-w 



(2) 



3+ ytJ 

r-,y,r- 

U.yJ 
3-,y,3L,^/ 



(28) 



and find 



= iET-[(r:t:^.-*^*i)-[-^* + ^^H^^ 

cr 

The solution will be sought by means of the functional form 

r[vl/t, M/] = f [M/t, + I E • [-9t + ^a{-idx)\ ■ 

for which the evolution equation reads as 



(29) 



(30) 



(31) 
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V. ANSATZ 

The general solution of the functional differential equation (31) is beyond our capacities and we have arrived at the 
only approximation one has to rely upon in this scheme, the truncation of the evolution equation. By assuming an 
analytical dependence in the field variables we write r[\l/^^, in the form of a functional Taylor series where arbitrary 
high powers of the fields appear and the coefficient functions, the IPI vertices, are arbitrary functions of their space- 
time or energy-momentum variables. The n-th order terms in this series represent n-body correlations and can be 
shown to be suppressed in the framework of the Wilsonian the external space renormalization group method in the 
limit when the sharp momentum space cutoff used approaches zero [11]. Though this result seems to be somehow 
fragile, the nonlocal interactions generated by the sharp cutoff imposed in the IR regime may substantially change 
the kinetic energy and thus invalidate the argument, we shall ignore the the terms n > 6 for the sake of simplicity. 
But the collinear divergences will bo taken more seriously and an arbitrary space-time dependence will be retained 
for the vertex functions. Our choice for the "potential" r[\E''l', is 

CT •^Xi,yi,X2,y2 

- / V f^ft ^\ ^j.„ + * .„ ^ ■„ ~) 

/ 'xi,yi,X2,J/2V*+;a:i ^+;X2 +;y2 ^ ^ -;xi ^ -;3/i ^ -;X2 .!/2 

Jxi,J/i,X2,J/2 



(32) 



where the effective vertices obey the symmetries ll^,y^,x2,v:, = 7x2,!/2,xi,yi = ~ixi,y2.x2,yi- 

We shall be interested in the question if the theory is strongly coupled aroiuid the Fermi points. It should be kept 
in mind that the eventual divergence of the vertex functions at some points does not necessarily mean that the model 
is strongly coupled in that kinematical regime. In fact, any physically acceptable definition of the strength of the two- 
body interaction involves a finite resolution in the energy-momentum space and seeks an average, a representative 
value of the four-point vertex function in a small but finite area. Such a construction can be approximated by 
the introduction of an appropriately chosen normalized, smooth smearing function Pp,q,r,s in the definition of the 
dimensionless running coupling constant. 
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^p—q-\-r—s,OPp,q,r,s'yp,q,r,i 



(33) 



where %,q,r,s is obtained from the vertex function by separating off the trivial energy-momentum conservation, 

/ Jpxi-igyi+irx2-isy2: - /o4^ 



7xi,!/l,X2,!/2 



Guided by such a construction we shall call the theory strongly coupled at a certain kinematical point if the four point 
function 7p,g,r,.s cither takes a imiformly large value around or displays a non-integrable singularity at that point. 

One can explore the weakly coupled regime of the theory in a consistent manner with our truncation, (32). But 
whenever the dimensionless running coupling constants assume values which is not anymore small with respect to one 
then our computation must stop. 



VI. EVOLUTION EQUATION 

In order to arrive at the evolution equation for our ansatz we need the inverse of the second functional derivative, 

/p(2) p(2) j.(2) j.(2) ^ 

\T/T >T/ . \T/T iT/ iT/T iT/T .T/T >T/T 



r(2) 



p(2) 

p(2) p(2) 
p(2) p(2) 



p(2) p(2) 

— jX ' +;y _ ^ 1 



.(2) 



.(2) 



.(2) 



^1 



.(2) 
.(2) 



(35) 



,Xj ^.y ,Xj — ;y/ 

computed in a manner which is consistent with the truncation. For this end we write it in the form 

r(2) =G'o^ + A'' +A" 



(36) 
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by means of the diagonal block matrix 



/G 


-1 




































\ 














composed of the free propagator 



the self energies 



I Z\,Z2,X,y^ —■,Zl^ — ,Z2 IZl,y,X,Z2 —;Zl^+,Z2 " /j;,2i,X,Z2 +1^2 \ 

-^x vl/t xEr vl/t \I/ , VEr vEr , f) 

/X,Zl,Z2,!/ +;Z2 'X,y,Zl,Z2 +;Zl^ +\Z2 lx,Zl,y,Z2 " 

" lZi,3/,Z2,a: — ;zi +;Z2 'zi,Z2,y,2; * — ;2;i ^ — ,Z2 ly,zi,Z2,x ^ +;Z2 ,^1 

" izi,x,y,Z2 —;zi +jZ2 iy,x,zi,Z2 +;zi +^^2/ 



\ lzi,x,Z2,y —■,zi''+;z2 



(37) 



(38) 



(39) 



and 



^x,y I 

J Z\^Z 



l'^ll,V,z„Z2'^+;z,-^+;z 







V 



-47;_,^,^_,^*L.,^*_;,j 



27;,zi,c.,z2*+;^i*+;^2 







2V' *t 

/zi,y,Z2,x +;2i ^ +;Z2 



2V' ^f'l' ^rt 

izi,y,Z2,x ^ — ;zi — ;Z2 



(40) 

The matrix elements of the second functional derivative in the desired approximation are finally obtained by means 
of the Neumann expansion, 

(r(2))-i = Go - Go • (A^ + A") • Go + Go • (A^ + A") • Go • (A>^ + A") • Go + • • • (41) 

which yields the evolution equation 

dxf = -^Tr{[-dt + e^-idx)] ■ [Go - Go • (A>< + A") • Go + Go • (A>< + A") • Go • (A>< + A") • Go]} (42) 

for our ansatz. The first contribution arising from Go in the second parenthesis on the right hand side is field 
independent and can be ignored. The quadratic and quartic terms generate the evolution for the self energy and the 
effective coupling constants, respectively, to be considered below. 



A. Quadratic part 

The quadratic order of the left and the right hand sides of the evolution equation are 

,Zl,Z2 -,Zl,Z2 ^ -,Z2/ 

J 2l,22 

and 



-T2 / [G+;z,,zA^ll + Al^)zuZ2G+-z2,zA-dt + e+i-idx)]z,,z, 

Jzi,Z2,Z3,Zi 

-T2 I [G-,z„zA^l2 + ^\2)z^,Z2G-,z2,zA-dt + e-{-idx)]z3,z, 

^ J Zi,Z2,Z3,Zi 

respectively. By equating these expressions one finds the evolution equations for the self energy, 

9aS-,x,j/ = / {G+;z4.,zilx,y,z^,Z2G+;z2,Z3[-dt + ^+{-idx)]z3,Zi 

^ J Zl,Z2,Z3,Zi 

+AG-.^^^^^i^^^^^^^^yG-.^^^^A-dt + e-{-idx)]z3,z^], 

d\^+,x,y = -T2 / {G-,Zi,zilx,y,zi,Z2G-,Z2,Z3[-(^t + (^-{-idx^s^Zi 

J Zi,Z2,Z3,Zi 

+4G+;2^,^i7;:^_^^_^_j^G+;22,z3[-^t + e+{-idx)]zs^^^}- 



(43) 



(44) 



(45) 



11 



The structure of these equations reflects the Callan-Symanzik strategy in the framework of evolution in A. In fact, 
the relation 

= = -^G. ■ l-a + ■ (46) 

where d'^ acts on the explicit, driving A-dependence produces 

d\^--X,y = d'x (7x,3/,Zl,22^+;Z2,Zl + 472i,Z2,X,3/<^-;Z2,Z3)) 



(47) 



J zx,zi 

the perturbative, one-loop Callan-Symanzik beta function for the propagator. The Fourier transform of the self energy 

S<.;.,,= / e^^f+^^«<5j,+,,oS.;p,„ (48) 
allows us to write the evolution equation in a simpler form, 

n V — f>' f ^'P,l':<i:<i ^ ^ 'I'-l'-'i-'l /AQ\ 

Jg j[~"t + ^t{'^"x)\ + 

It is worthwhile noting the appearance of different derivatives on the two sides of these equations. In fact, had we found 
d\ on both sides the solution would have produced the standard one-loop results for the self energy. The derivative 
d'^ on the right hand side makes sure that the effective vertices 7^ , 7" which in principle depend on A are nevertheless 
kept intact. This procedure shows the essence of the infinitesimal blocking version of the renormalization group 
strategy: the for an infinitesimal change of the control parameter can safely be computed perturbatively by means of 
the actual effective vertices and the nonperturbative resummation arises from the repetition, i.e. the integration of 
this procedure. 



B. Quartic part 

The evolution equation gives in the quartic order 

A 5a 7^1 ,yi,X2,y2 ~ / {G+-ze,Zi{ ^lxi,yi,Zi,Z2^+;Z2,Z3^Z3,Z4,,X2,y2 '^^Zi,Z2,X2,y2^+'^-2,Z3lxi,yi,Z3,Z4 

J Zi,...,Zq 

'^lxi,Z2,Zi,y2^-;Z2,Z3lz3,yi,X2,Zi ~^ lz2,yi,Zi,y2^-\Z3,Z2lxi,Z3,X2,Zi)^+;Zi,Z5 [~dt + e-|-(~i5x)]z5,25 
~^^—'jZ6,Zi ilzi,yi,X2,Z2^+'tZ2,Z3^xi,Z4,Z3,y2 ~^ Tzj ,yi ,Z2 ,1/2 ^+'^3 Txi ,24 ,a;2 ,Z3 

~^^Zi,Z2,Xi,yi^-;Z2,Z3l'z3,Z4,X2,y2 ~ ^1' Zi ,Z2 ,X2 ,y2^ -y^2 ,Z3l Z3,Z4,Xl ,yi)^ -■,Z4,Z!i [~ + £- { — idx)]z!i ,Ze} 

,yi,X2,V2 ~ / {^+;Z6,Zl{ 'yxi,yi,Zl,Z2^+'>Z2,Z37x2,V2,Z3,Z4 7x2 , J/2 ,21 ,22 ^+;22 ,237x1 ,1/1 ,23 ,24 

J Zl,...,Z6 

'^^Xi,y2,Zi,Z2^+\Z2,Z3lx2,yi,Z3,Z4 + 7a;2 , J/1 ,2i ,22 ^-|-;22 ,23 7xi , J/2 ,23 ,24 ) ^-|-;Z4 ,25 [~dt + C-F (" ^^^X )] 25 ,26 
"1"^— ;26,Zl {~^^Jzi,Z2,Xi,yi ^—;Z2,Z3l Z3,Z4,X2,y2 ~ Zi ,Z2 ,X2 ,y2^ ,Z3^ Z3,Z4,Xl ,yi 

~^^^^Zl,Z2,Xl,y2^—lZ2,Z3^Z3,Z4,X2,yi ~^ 1^721 ,22 ,X2 ,yi ;Z2 ,23 723 ,24 ,Xl , J/2 (^0) 

~4')'z2,J/l,2l,J/2^-;2:3,Z2 7xi,Z4,rr2,Z3 + ^lz2,y2.Zi,yi^-;Z3,Z2lxi,Z4,X2,Z3)^Z4,Z!j [~dt + C- (" ^t^X )] 26 ,26 } ' 

The Callan-Symanzik strategy can be used again to simplify the right hand sides as 

,yi,x2,y2 A / L"^-l-;24,2i 1 ,yi ,21 ,22 ^+;z2 ,23723 ,24 ,X2 ,y2 ,yi,2l,y2^~;^3,22 7xi,23,X2,24) 

J Z\^. . .,24 

"1"^— ;24,2l (72l,yi,X2,22^-|-;22,23 7xi,24,23,J/2 ~ ^721 ,22 ,X2 ,J/2 ^^2 ,23 723 ,24 ,Xl ,yi )] 

5a7xi ,yi,X2,y2 ~ j [G+-^Z4,Zi{ 7x2, J/2, 21, 22^+;2:2, 23 7x1, yi, 23, 24 + 7x2,2/1, 21, 22^+:2:2, 23 

7xi,y2,23,24) 

J 21 ,...,24 

;24,2l (~1672l,22,Xi,2/i ;22,23 723,24,X2,J/2 "'"^721 ,22 ,Xl ,3/2 ;«2 ,23723 ,24 ,X2 ,yi 

-4722 ,!/l,^l,!/2^-;Z3,22 7xi,24,X2,23 )]• (51) 
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where the one-loop beta functions appear on the right hand side. These relations assume the form 

^\lp,q,r,s ~ I [^+;M(~7r,s,u,u^+;u7p,q,v,M + 7r,g,u,u^+;u7p,s,v,M) 

'^(^—;u{~^Q^u,v,p,q^—;v^v,u,r,s ^^lfu,v,p,s^ —;v^v,u,r,q ~ '^7ti,q,u,s^— ;w7p,M,r,u)]' (^^) 

in momentum space. 

VII. RUNNING COUPLING CONSTANT AT THE FERMI POINT 

We first confront our evolution equation with the standard one-loop Wilsonian renormalization group approach 
to the Luttinger model [12, 13] where the effective bare coupling constant is supposed to be energy-momentum 
independent and its value is given by the coefficient of the quartic term in the formion fields taken at the one-shell 
Fermi-point. Furthermore we set 7" = at the initial conditions and at any subsequent stage of the evolution. The 
evolution of the self-energy, Eq. (49) is now 

r -""(Air"'''" 

Because the sclf-cncrgy corrections are 0{"f), S± can safely be ignored in the right hand side in the asymptotically 
weakly coupled initial conditions, imposed at Aq « 0. The right hand side starts with a p-independent value at the 
beginning and obviously keeps this simple form during the rest of the evolution. 



a.E,., = -^„x(A)/^-,-— ^. (54) 



g 



The self energy remains real and it can be absorbed in the renormalization of the chemical potential. 
Let us continue with the evolution of 'y^{X) which, according to Eq. (52) is governed by the equation 



dxrW=d'xf 7^'(A)[G + ;„G_;_„ + (A)G_;„G+;„ + 2fc,]. (55) 



The more detailed form of the two contributions on the right hand side, 

rA r+00 ^^0 1 



, 2n {ivP - u) {ivP + u) 

J - 2^y_ -2^WTu)W^y ^^^^ 

shows the vanishing of the beta function, 

dxl^ = 0, (57) 

the absence of dressing in agreement in Refs. [12, 13]. 

Notice that the crucial step of this argument, the formal cancellation of the two integrals in Eqs. (56), takes place 
at the Fermi points only. A more careful study of the evolution equation (52) presented in the next Section reveals 
that the incomplete cancellation between the contributions on the right hand side of Eq. (55) generates strongly 
coupled dynamics in the vicinity of the Fermi points. 



VIII. 15 FUNCTIONS OFF THE FERMI POINT 



Let us consider now consider the complete evolution. One side of the problem is the evolution of the propagator, 
driven by the quadratic part, Eq. (49). The other one is the evolution of the two-body interactions, given by the 
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quartic pieces, Eqs. (52). This solution of this set of coupled integro-differential equations represents an involved 
numerical problem, owing to the fact that the evolution of momentum dependent functions actually involves infinitely 
many parameters. We reduce the complexity in this work in order to arrive at an analytically treatable problem by 
suppressing the evolution of the quadratic part in our ansatz (32), i.e. we impose S = 0, and leave the building up of 
the true interactions for the quartic effective vertices. 

The quartic evolution, shown by Eqs. (52), is still a formidable numerical problem due to the singular structures 
generated in the integrand and will be simplified in the following manner. We assume that the naive ideas aboTit the 
dominance of the interaction vertices by momentum independent parts is valid and introduce the usual running cou- 
pling constants, given by the strength of effective interactions at a certain kinematical point. The only generalization 
with respect to the previous Section is that the running coupling constant can be defined off the Fermi points. We 
shall then check the consistency of this picture by looking into the magnitude and the momentum-dependence of the 
effective two-body vertices obtained in this approximation. 

The /3-functions for the coupling constants J^^g^f^g = 7^ and j'^^g^f^g ~ 7" ^ given point (p, q, f, s) are introduced 

as 



f3" 



(58) 



Let us assume that the region around the point p, q, r, s dominates the dynamics. Then we can approximate the solu- 
tion of the integro-differential equations (52) by replacing the vertex function 7p,q,r,s by the momentum-independent 
7 = 7p,Q,f ,s in the loop integrals. The approximation is better if the point {p, q, f, s) is closer to the Fermi point. 
We shall parametrize the energy-momentum as 



(-iw, p) = i±p + Ap,p±kF) 



(59) 



where the plus and minus sign stand for the right and left moving particles, respectively. The parameters Ap, Ag, 
Ar, and measure the distance from the mass shell. The momentum conservation yields s — r = p — q and the 



energy conservation leads to the condition 2(p — q) = Ap ■ 
for 7^ and 7", respectively. 

The beta functions in Eqs. (52) can then be written as 



Aq + Ar 



A,, and 



Ag + Ar 



Ag in equations 



/?" = -A2(/3,fx^x2 + /3;;"Y.2^, 



(60) 



with 



-27rln 



{2q -Ag + Ar){2f -Ag + Ar)[{2r + Ap + Arf - 4A^] 



{2r + Ap + Ar){2p -Ap- A,)[(2q + A, - A^Y - 4A2] 
p-q p-q 




Ag + 2q-2kF , 
r — q + 2kF 



p + r + 2kF 



Ag-Ap + 2{q - p) 



- Aq + 2q- 2kF Ap + Ar + 2r + 2kF 



(61) 



Sec Appendix B for details. 

The beta functions display singularities arbitrarily close to the Fermi points which lead to singularities in the 
effective coupling constants [29]. We shall integrate out the beta functions in order to see the singularity structure 

in a bit more detailed manner. Two kinematical regions will be considered, (a): q = p, A = 0, f,p <C ^FjA, and 



(b): q = p, 2f + Ap + Ar = 2p + Ag — A^, r,p <^ kF, where /3J 
0{r/kF) direction dependent terms and 



0, a: 



47r, /3|'|" ~ 367r by ignoring 0{p/kF), 



X X 
X 



2tt^ 



A2 



27rln 



2p-Ag+Ar 



(6). 



(62) 
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The evolution of 7^ is autonomous 

dX 



^5^ = -^'/3^'7"^ (63) 



with the solution, 



where 7^(0) = 7(^ denotes the bare coupling constant which is supposed to be small. |7(^| <C 1. The solution (64) 
shows that our truncation resunis the geometrical series of the diagonal part of the vertex function. The truncation 
of the effective action is consistent and the system remains weakly coupled for (i^^ > 0, in which case |7^(A)| 
decreases during the evolution. On the contrary, for 7^ < — 1 a Landau-pole is predicted at 



before the desired bare coupling strength 70 is reached. The Landau-pole appears in region (a) due to the cutoff 
effects when the inequality 

is satisfied by the bare theory. This pole would be absent in a typical high energy application when the cutoff is 
large but may appear as a non-universal effect when the model is considered as an effective theory with low cutoff, 
A <C fci?. The Landau-pole is always present in region (b), in the vicinity of the mass shell and the Fermi point, where 
the inequality 



2p 



(67) 



can always be satisfied for an arbitrary value of 7,^ . 

The effective coupling constants 7" is obtained by integrating the beta function 

/3" = -4A27r(7><2+97"2). (68) 
The evolution is autonomous for 7^ (A) ~ and gives 



with another Landau-pole when 



^0 < 



When 7^ (A) approaches a Landau-pole then 7" —00 because the singularity is non-integrable. 

It is worthwhile noting that fi^ ^ and together with it the coupling strengths 7^ (A) and 7"(A) are scale invariant in 
regime (b), i.e. the external scale parameters kp and A decouple from the dynamics. As a result, the beta functions, 
obtained by means of the usual renormalization group schemes where the energy- momentum components of subtraction 
point scale with A in a homogeneous manner, are vanishing. Furthermore, the beta functions (61) and the effective 
coupling strengths become scale invariant in the limit A ^ 00 at the subtraction point, ujp = ujg = uij- = ojg = 0, 
p/3 = q = -f = 2, used in Ref. [10]. 



IX. SUMMARY 



The transmutation of quasiparticles between the short and the long distance scales is a challenging problem from 
particle to condensed matter physics. We considered one dimensional massless fermions, the Luttinger model in this 
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work and applied the functional form of the internal space renormalization group with sharp momentum cutoff which 
yields a functional differential equation for the effective action. The initial conditions for the differential equation 
are known. This equation was truncated in such a manner that the weakly coupled nature of the interactions can be 
checked. The solution shows strong interactions and singular momentum dependence in the vicinity of the Fermi points 
in the two-body channels and Landau poles for arbitrary bare parameters. But our results are rather preliminary, one 
should systematically enlarge the ansatz for the effective action and study the stability of Landau poles in the result. 

We argued that cutoff effects, in particular Landau poles, are generated by the C(A~^) part of the renormalized 
action in the normal ordered schemes, such as bosonization. We believe that more work is needed on the one hand, to 
increase the reliability of the results obtained by completely regulated bare theories and on the other hand, to clarify 
the role the 0{A~'^) interactions may play in the bosonized theory. 
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APPENDIX A: NOTATION, CONVENTIONS 



The space-time and Fourier-space integrals are defined as 



(Al) 



where x = {t,x), p = {(jJ,p). In Euclidean space-time one has px = p'^x'^ + px. The letters a, b, . . . and p, q, ■ ■ ■ 
denote space-time or energy-momentum variables, respectively. The integration is sometime shown as scalar product, 
f -9 = J^fxgx- 

The functional derivative is defined in d-dimensional space-time lattice as 



6^^ a^at 

where ^l/-^ is the lattice field variable whose dimension is removed by the lattice spacing and the factor l/a^at is 
needed in order to satisfy the relation 

t|- / /(<^x) = /'(*,). (A3) 

In an analogous manner we have 



APPENDIX B: ONE-LOOP INTEGRALS 



This Appendix contains some details concerning the computation of the one-loop integrals of Eqs. (52) by means 
of sharp momentum space cutoff. In order to simplify the results we assume that the absolute value of the momenta, 
counted from the Fermi point does not exceed A/4. 



For the first equation we need tlie following integrals, 

^3 = ^ J G--tG+-t-q+r, 



G--tG--t+i — s) 



where 



Go 



A 



-dt + e„{-idx) ' 

One uses the residuum theorem in integrating over the frequency in the integral 



/oo r+oo 
du / 
-OO J — 



h 



du!- 



OO — oo 



and finds 



which results 



/i = 2i7r / du 



{lo + iu) {uj + Up - Uq + i{u + p - q)) 

Pa{u)pa{u+p - q)[(d{-u) - Q{-u + q-p) 
oJp -uJq +i{p- q) 



J ^ 27r(p - q) 
^ Ap - Ag + 2{p - q) ■ 

The other three integrals can be computed in a similar manner. For 

''^'^ pA{u)pA{u-p-r) 



h = 



j du i 

J —CO J —CO 



du)- 



we find 



for \p\, \r\ < A/4 where 



{lj + iu) [w — Wp — LOr — i{u — p — r)) 

h = ^2(0) - F2{-K) - [F2{K) -F2{p + r)] 



F2{u) = — Trln 



P + r .LOp + UJr 

u t—^ 



-In 



u — r 



Ap + Ar 



The final result is 



-In 



(r+^){p-^) 



The integral 



/oo f+CO 
du / du} 
-co J — CO 



Pa{u)pa{u -q + r) 



{oj — iu) {oj — (Jq + (Jr + i{u — q + T — 2fcir)) 



becomes 



/3 = -F3(A) - Fsi-A) + F3{0) +F3{q- f) 
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where 

^3(11) = Trln 

after the integration over the frequency. Trivial steps give 



, r-q .ujg~ujr 



FinaUy, the last integral, 



gives 



I3 = TT In 



00 r+oo 



I du I 

J — 00 J —00 



du: 



Pa{u)pa{u - s + r) 



{uj — iu) {u + LOr — LOs + i{—u — r + s)) 



7. 



2i7r / du 

-A 



e(n) 



+ ■ 



@{u + r — s) 
— + i{s — r) ' — — i{s — r) 
27r(s — f) 



As-Ar + 2{s - r) ■ 



(B12) 



(B13) 



(B14) 



(B15) 



2. 7" 



For the second equation in Eqs. (52) we need the following five integrals, 

•^^ — 1^ J G+-tG+-r-s+t, 



J2 

J3 



— ^ J G+-tG+-r-q+t, 



1 



G--^tG--p+r-t, 



"^4 — / G--^tG--^t-r+si 



G-,tG 



(BI6) 



Due to the left and right symmetry properties J4 = Ji and J5 = J2. The same procedure than in the case of the 
previous integrals gives 



Ji 



/oo p-\-oo 
du / ( 
-00 -/ —00 



Pa{u)pa{u- s + r) 



(w + iu) {uj + LOr — LOs + i{u + r — s)) 



2i'ir / dupA{u)pA{u — s + r) 

—00 

27r(f — s) 



6(— u + s — r) 

— Ws + — s) ' LOs — UJr — ^(^ — s) 



e{-u) 



+ 



and 



J2 



- A, + 2(f - s) ' 



(B17) 



Pa{u)pa{u -q + r) 



{lo + iu) {lo + ujr — ijJq + i{u + T — q)) 



= 2n 



2in / dupA{u)pA{u — s + r) 

J —00 

r — q + 2kp 
Ar - A„ + 2q- 2kF 



e{-u) 



+ 



Q{—u — r + q) 



<r - cjq + i{r - q) ujq-ujr-i{r-q) 



(B18) 
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Finally, the last independent integral is 



Js = - 
= -2i7r 
= -27r 



du 



Pk{u)pk{u -p-r) 



{lj — iu) {u) — uip — uir — i{u — p~ r — 2kF)) 

e(M) 



dupK{u)pA{u -p-r) 

OO 

p + r + 2ki.- 



Q{u — p — r — 2k[ 



i{—p — r — 2kp) Up + ujr + i{—p — r — 2kp) _ 



Ap + Ar + 2f + 2kF 



(B19) 
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